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In this paper, we compare some popular dark energy models under the assumption of a flat universe by using 
the latest observational data including the type la supernovae Constitution compilation, the baryon acoustic 
oscillation measurement from the Sloan Digital Sky Survey, the cosmic microwave background measurement 
given by the seven-year Wilkinson Microwave Anisotropy Probe observations and the determination of Hq from 
the Hubble Space Telescope. Model comparison statistics such as the Bayesian and Akaike information criteria 
are applied to assess the worth of the models. These statistics favor models that give a good fit with fewer 
parameters. Based on this analysis, we find that the simplest cosmological constant model that has only one 
free parameter is still prefened by the current data. For other dynamical dark energy models, we find that 
some of them, such as the a dark energy, constant w, generalized Chaplygin gas, Chevalliear-Polarski-Linder 
parametrization, and holographic dark energy models, can provide good fits to the current data, and three of 
them, namely, the Ricci dark energy, agegraphic dark energy, and Dvali-Gabadadze-Porrati models, are clearly 
disfavored by the data. 
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I. INTRODUCTION 

Dark energy has become one of the most important issues of the modern cosmology ever since the observations of type la 
supernovae (SNe la) first indicated that the universe is undergoing an accelerated expansion at the present stage [1]. However, 
hitherto, we still know little about dark energy. The limited information we know about dark energy includes: it causes the 
cosmic acceleration; it accounts for two-thirds of the cosmic energy density; it is gravitationally repulsive; it does not appear 
to cluster in galaxies; and so on. Many cosmologists suspect that the identity of dark energy is the cosmological constant that 
fits the observational data well. While, one also has reason to dislike the cosmological constant since it always suffers from 
the theoretical problems such as the "fine-tuning" and "cosmic coincidence" puzzles [2]. The fine-tuning problem, also known 
as the "old cosmological constant problem," is motivated by the enormous discrepancy between the theoretical prediction for 
the cosmological constant and its measured value. The so-called "new cosmological constant problem," namely, the cosmic 
coincidence problem, questions why we just live in an era when the densities of dark energy and matter are almost equal, which 
also indicates that the cosmological constant scenario may be incomplete. Thus, a variety of proposals for dark energy have 
emerged. 

The possibility that dark energy is dynamical, for example, in a form of some light scalar field [3], has been explored by 
cosmologists for a long time. A basic way to explore such a dynamical dark energy model in light of observational data is to 
parameterize dark energy by an equation-of-state parameter w, relating the dark energy pressure p to its density p via p - wp. 
In general, this parameter w is time variable. The most commonly used forms of w(a) involve the constant equation of state, 
w = const., and the Chevalliear-Polarski-Linder form [4], w{a) - wo + (1 - a)wa, where wo and Wa parameterize the present-day 
value of w and the first derivative. There are also many other dynamical dark energy models which stem from different aspects of 
new physics. For example, the "holographic dark energy" models [5-1 1 ] arise from the holographic principle of quantum gravity 
theory, and the Chaplygin gas models [12-14] are motivated by brane world scenarios and may be able to unify dark matter and 
dark energy. In addition, there is also significant interest in modifications to general relativity, in the context of explaining the 
acceleration of the universe. The Dvali-Gabadadze-Porrati models [15-17] arise from a class of brane-related theories in which 
gravity leaks out into the bulk at large distances, leading to the accelerated expansion of the universe. 

In the face of so many competing dark energy candidates, it is important to find an effective way to decide which one is 
right, or at least, which one is most favored by the observational data. Although the accumulation of the current observational 
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data has opened a robust window for constraining the parameter space of dark energy models, the model filtration is still a 
difficult mission owing to the accuracy of current data as well as the complication caused by different parameter numbers of 
various dark energy models. In this paper, we make an effort to assess some popular dark energy models in light of the latest 
observational data, including the Constitution SN data [18] and other cosmological probes such as the distance information 
measured by Wilkinson Microwave Anisotropy Probe (WMAP) [19], and the Baryon Acoustic Oscillations (BAO) [20, 21]. To 
make a comparison for various dark energy models with different numbers of parameters and decide on the model preferred by 
the current data, following Ref. [22], we apply model comparison statistics such as the Bayesian information criterion (BIC) 
[23] and the Akaike information criterion (AIC) [24] in our analysis. 

This paper is organized as follows. In Sec. II, we discuss the information criteria in the context of dark energy model selection. 
In Sec. Ill, we give details of the observational data sets used. Section IV describes nine popular dark energy models and assesses 
which one is preferred by the current data. The results are discussed in Sec. V. 



II. METHODOLOGY 

In this work we employ the;if^ statistics. For a physical quantity ^ with experimentally measured value ^obs, standard deviation 
CTf , and theoretically predicted value f,;,, thex^ value is given by 

• (1) 



The total is the sum of all;\^^s, i.e.. 



- (2) 



The observational data we use in this paper include the Constitution SN la sample, the Cosmic Microwave Background (CMB) 
measurement given by the seven-year WMAP observations, the BAO measurement from the Sloan Digital Sky Survey (SDSS), 
and the measurement of Hq from the Hubble Space Telescope (HST) [25]. 

However, the statistic alone cannot provide effective way to make a comparison between competing models since this 
method is based on the assumption that the underlying model is the correct one. The statistics is good at finding the best-fit 
values of parameters but is insufficient for deciding whether the model itself is the best one. Since in general a model with 
more parameters tends to give a x^,,i„, it is unwise to compare different models by simply considering with likelihood 
contours or best-fit parameters. Instead, one may employ the information criteria (IC) to assess different models, which is also 
based on a likelihood method. In this paper, we use the BIC [23] and AIC [24] as model selection criteria. According to these 
criteria, models that give a good fit with fewer parameters will be more favored. So, these criteria embody the principle of 
Occam's razor, "entities must not be multiplied beyond necessity." The applications of the BIC and AIC in a cosmological 
context can be found in, e.g., Refs. [26, 27]. 

The BIC, also known as the Schwarz information criterion [23], is given by 

BIC = -2 In + k\nN, (3) 

where Xmn.v is the maximum likelihood, k is the number of parameters, and is the number of data points used in the fit. Note 
that for Gaussian en:oK,x\i,j - -2 InXmat, and the difference in BIC can be simplified to ABIC - ^xlu,, + ^klnN. A difference 
in BIC (ABIC) of 2 is considered positive evidence against the model with the higher BIC, while a ABIC of 6 is considered 
strong evidence. 

The AIC [24] is defined as 

AIC = -21ni:,„„ + 2^, (4) 

which gives results similar to the BIC approach, but it should be pointed out that the AIC is more lenient than BIC on models 
with extra parameters for any likely data set In > 2. Also, in this case, the absolute value of the criterion is not of interest, only 
the relative value between different models, AAIC = ^Xmm + 2AA:, is useful. As mentioned in Ref. [28], there is a version of the 
AIC corrected for small sample sizes, AIC^ = AIC + 2k{k - l)/iN - k - 1), which is important for N/k < 40. Obviously, in our 
case, this correction is negligible. 

It should be noted that the information criteria alone can at most say that a more complex model is not necessary to explain 
the current data, since a poor information criterion result might arise from the fact that the current data are too poor to constrain 
the extra parameters in this complex model, and it might become preferred with improved data. Actually, this is just the current 
situation for dynamical dark energy models. 
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A more sophisticated method for model selection is provided by the Bayesian evidence which does not simply count param- 
eters, but considers how much the allowed volume in data space increases due to the addition of extra parameters, as well as 
any correlations between the parameters. So, the Bayesian evidence requires an integral of the likelihood over the whole model 
parameter space, which may be lengthy to calculate, but avoids the approximations used in the information criteria and also 
permits the use of prior information if required. This method has been applied in a variety of cosmological contexts; see, e.g., 
Refs. [29-31]. Information criteria require no assumptions for the prior or the metric on the space of model parameters. In this 
paper, we will use the first approximation provided by the information criteria without calculating the full Bayesian evidence. 
This simpler version is sufficient for our purpose. 



III. CURRENT OBSERVATIONAL DATA 



In order to test the different dark energy models, we have used the observational data currently available. In this section, we 
describe how we use these data. 



A. Type la supernovae 

Up to now, SNe la provide the most direct indication of the accelerated expansion of the universe. It is commonly believed 
that these SNe la all have the same intrinsic luminosity, and thus they are used as "standard candles." Therefore, measuring both 
their redshift and their apparent peak flux gives a direct measurement of their luminosity distance di as a function of redshift 
z. The function di{z) encodes the expansion history of the universe so that by which the information of dark energy can be 
extracted. 

In this paper, for the SN la data, we use the Constitution sample including 397 data that are given in terms of the distance 
modulus fiobsizi) compiled in Table 1 of Ref. [18]. The theoretical distance modulus is defined as 

m(zi) = 5 logio Diizi) + fio, (5) 

where /yq = 42.38 - 5 logjy /; with /? the Hubble constant //q in units of 100 km/s/Mpc, and the Hubble-free luminosity distance 
Dl = Hodi is 

Ddz) = ^7^sinn( f (6) 

vin^ V Jo Eiz')/ 

where E{z) = H{z)/H{), is the fractional curvature density at z = 0, and 

Sinn ( y/l^hi'c) ( ViH^x)/ VjO^, if < (or ^ = +1), 
^ ' = \ X, if = (or ;t = 0), 



sinh( y/\n[\x)/ Vl?¥, if £2^ > (or k ^ -1). 



Thex^ for the SN data is 



where fiobsizi) and cr, are the observed value and the corresponding 1 cr error of distance modulus for each supernova, respectively, 
and 9 denotes the model parameters. 

What should be mentioned is that since the absolute magnitude of a supernova is unknown, the degeneracy between the Hubble 
constant and the absolute magnitude implies that one cannot quote constraints on either one. Thus the nuisance parameter Hq 
in SN data is not the observed Hubble constant and is different from that in the BAO and CMB data. Therefore we should 
analytically marginalize over Hq in the SN data. 

Following Refs. [32], the minimization with respect to yUo can be made trivially by expanding the;^'^ of Eq. (7) with respect to 
yuo as 

xl^iO) = A(e) - + f^lC, (8) 



where 

[i^ohs(zi) - i^ihizr, Mo = 0, e)] 

cr2 



^^g-j ^ y inohsizi) - fiihizr, fio ^ 0,0)] 
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err 

Evidently, Eq. (7) has a minimum for jiQ = B/C at 

9 B(0)- 

xl^(e)^A(e)-^. (12) 

Since = XsNmin^ instead minimizing we will minimize j^^^ which is independent of the nuisance parameter fiQ. 



B. Baryon Acoustic Oscillations 



Aside from the SN data, the other external astrophysical results that we shall use in this paper for the joint cosmological 
analysis are the BAO and CMB data. First, we describe how we use the BAO data. The BAO is a powerful probe of dark energy, 
since it can be used to measure not only the angular diameter distance Da(z) through the clustering perpendicular to the line of 
sight, but also the expansion rate of the universe H{z) through the clustering along the line of sight. However, the current data 
are not accurate enough to allow us to extract D^iz) and H{z) separately. Actually, the BAO currently can barely be measured in 
the spherically-averaged correlation function [33]. 

The spherical average gives us the following effective distance measure [34] 



Dviz) = 



il+z)-Dl(z) 



H(z) 



1/3 



where D^iz) is the proper (not comoving) angular diameter distance, 

1 sinn( ^Mk\Ho £ 4^)1 ^IWk 



Da{z) 



Ho 



a + z) 



(13) 



(14) 



The BAO data from the spectroscopic SDSS Data Release 7 (DR7) galaxy sample [21] give Dv(z = 0.35)/Dv(z = 0.2) 
1 .736 + 0.065. Thus, the for BAO data is. 



2 _ ( Dviz ^ 0.35)/Dv{z ^ 0.2) - 1.736 
^^^0 " I 0.065 



(15) 



C. Cosmic Microwave Bacliground 

The CMB is sensitive to the distance to the decoupling epoch via the locations of peaks and troughs of the acoustic oscillations. 
In this paper, we employ the "WMAP distance priors" given by the seven-year WMAP observations [19]. This includes the 
"acoustic scale" U, the "shift parameter" R, and the redshift of the decoupling epoch of photons z,- 

The acoustic scale 1^ describes the distance ratio D^izt)/ r^iZt), defined as 

,,.(1,,.)!^, (16) 

where a factor of (1 + Zt) arises because Z)a(z*) is the proper angular diameter distance, whereas rj(z,) is the comoving sound 
horizon at z,. The fitting formula of rs{z) is given by 



da_ 

V3 Jo a^H{a) ^1 +(3Qfe/40y)fl' 



rAz)^-\= f , „ „ , (17) 

V3 Jo 



where Q.b and are the present-day baryon and photon density parameters, respectively. In this paper, we fix Qy - 2.469 x 
lO^^h^^ (for Tc,„h = 2.725 K) and Q/j = 0.022765/1"^, which are the best-fit values given by the seven-year WMAP observations 
[19]. We use the fitting function of z» proposed by Hu and Sugiyama [36]: 

z, = 1048[1 +0.00124(Qi/iV™][l -t-g^n^/iV^], (18) 



5 



where 



0.0783(Qa/iO 



2\-0.238 



1 +39.5(fii,/l2)0.763' S2 1 +21.1(Qft/l2)l-81- 

The shift parameter R is responsible for the distance ratio DAiZt)/H^\zt), given by [37] 



0.560 



(19) 



(20) 



Actually, this quantity is different from Da(z*)/H^^(z*) by a factor of VI + z*, and also ignores the contributions from radiation, 
curvature, or dark energy to H{zt). Nevertheless, we still use R to follow the convention in the literature. 

Following Ref. [19], we use the prescription for using the WMAP distance priors. Thus, the;^'^ for the CMB data is 



maximum likelihood values; Ia{z») - 302.09, R{Zt) 
Ref. [19]: 



trix. The se 

1091.3. The inverse covariance matrix is also given in 



where jc, - {Ia,R, z*) is a vector, and (C '),; is the inverse covariance matrix. The seven-year WMAP observations [19] give the 



1.725, and: 



2.305 29.698 -1.333 ^ 
29.698 6825.27 -113.180 
-1.333 -113.180 3.414 



(22) 



D. Hubble constant 



In this paper we also use the prior on the present-day Hubble constant, Hq - 74.2 + 3.6km/s/Mpc [25]. In Ref. [25], the 
authors obtain this measurement result of Hq from the magnitude-redshift relation of 240 low-z type la supernovae at z < 0.1. 
The absolute magnitudes of these supernovae are calibrated by using new observations from HST of 240 Cepheid variables in 
six local type la supernovae host galaxies and the maser galaxy NGC 4258. It is remarkable that this Gaussian prior on //q has 
also been used in the analysis of WMAP 7-year observational data [ 19] . The function for the Hubble constant is 

2 / /I- 0.742 \- 



E. Combining the constraints 



Since the SN, BAO, CMB and Hq are effectively independent measurements, we can combine our results by simply adding 
together the;if2 functions. Thus, we have 

^xIn+ xIao + xImb + xl ■ (24) 
Note that jf^^ andx^gAO '^^^ ^^^^ "^^^^^ ^cmb ^"'^^a '^^^ ^'■^^^ relevant to h. 



IV. DARK ENERGY MODELS 



In the standard homogeneous and isotropic Friedmann-Robertson- Walker (FRW) universe, the Friedmann equation is ex- 
pressed as 

, , 3Mi,k 

mliH^ = p - — f-, (25) 

where Mpi si/ y/SnG is the reduced Planck mass, p is the total energy density containing contributions from cold dark matter, 
baryons, radiations, and dark energy, and k describes the spatial geometry of the universe. Usually, the properties of dark energy, 
such as its equation of state parameter w, are degenerate with the spatial curvature of the universe fl^. Thus, although in principle 
one should include as an additional parameter when fitting dark energy models in light of observational data, the current data 
are not accurate enough to distinguish between w{z) and Qt, owing to the degeneracy of them. On the other hand, it is well 
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known that most inflation models in which the inflationary periods last for much longer than 60 e-folds predict a spatially flat 
universe, Q^. ~ 10"^. Actually, the inflation theory has become a paradigm in the modern cosmology and it has received strong 
support from the CMB observations. Under such circumstances, therefore, in this paper we shall use a "strong inflation prior," 
imposing a flatness prior, and explore dark energy models in the context of such inflation models. 
In a spatially flat FRW universe (O^ = 0), the Friedmann equation (25) reduces to 

3Ml,H^ = p,„(l + z? + p,(l + zf + Pdem(z), (26) 

where p„,, pr and pdeiO) are the present-day densities of dust matter, radiation and dark energy, respectively, and f{z) = 
Pdeiz) / pde{0) is givcn by the specific dark energy models. This equation is usually rewritten as 

E{z) = H{z)/Hq = [q,„(1 + zf + a(l + zf + (1 - Q,„ - ^r)f{z)f'^ . (27) 

Note that the radiation density parameter Q,- is the sum of the photons and relativistic neutrinos [35], 

Q, = + 0.2271A^,//), (28) 

where N^ff is the efl^ective number of neutrino species, and in this paper we take its standard value 3.04 [19]. In some cases, the 
evolution of the dark energy density parameter HdAz) = Pde{z)/OMpiH^) is determined by a differential equation, and thus one 
should express the Friedmann equation as 

In what follows, we choose nine popular dark energy models and examine whether they are consistent with the data currently 
available to us. We divide these models into five classes: 

1. Cosmological constant model. 

2. Dark energy models with equation of state parameterized. 

3. Chaplygin gas models. 

4. Holographic dark energy models. 

5. Dvali-Gabadadze-Porrati (DGP) brane world and related models 

It should be mentioned that the dark energy models discussed in this paper involve those invoke a variation in the equations 
governing gravity, say, the DGP model. Actually, in some cases the two are interchangeable descriptions of a single theory, 
albeit they might affect the structure growth in different manners (see, e.g., Refs. [38, 39]). In this paper, we ignore the exiguous 
difference between them, and call them uniformly the "dark energy models." The models discussed and the parameters that 
describe each model are summarized in Table I. Note that when fitting the data there is an additional parameter, h, which is not 
considered as a model parameter but is included in k when calculating AIC and BIC. The fit and information criteria results are 
summarized in Table II. 

Next, we shall outline the basic equations describing the evolution of the cosmic expansion in each of the competing dark 
energy models, calculate the best-fit values of their parameters, and find their corresponding AAIC and ABIC values. 



A. Cosmological constant model 

The cosmological constant A was first introduced by Einstein [40] with a wrong motivation, but nowadays it has become the 
most promising candidate for dark energy responsible for the current cosmic acceleration. While it has been suffering from the 
theoretical problems, it can explain the observations well. The cosmological model containing a cosmological constant and cold 
dark matter (CDM) component is usually called ACDM model in the literature. The unique feature of the cosmological constant 
is that its equation-of-state parameter w has the value -1 at all times, so in this model we have 

E{Z) = V"m(l + zf + a(l + Z)'* + (1 - - a). (30) 

It is obvious that this model is a one-parameter model, with the sole independent parameter Q,„ . 

The best-fit values of parameters (including the model parameter Q„, and the dimensionless Hubble parameter h) and the 
corresponding are: 



Q,„= 0.275, /i = 0.704, 468.461. 



(31) 
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TABLE I: Summary of models 



Model 


Abbreviation" 


Model parameters" 6 


Number of model parameters (/: - 1) 


Cosmological constant 


A 




1 


Constant w 


w 


n„,, w 


z 


Chevallier-Polarski-Linder 


CPL 




3 


Generalized Chaplygin gas 


GCG 


As, a 


2 


Holographic dark energy 


HDE 


^m, C 


2 


Agegraphic dark energy 


ADE 


n 


1 


Ricci dark energy 


RDE 


n,„, a 


2 


Dvali-Gabadadze-PoiTati 


DGP 


Clfii 


1 


Phenomenological extension of DO? 


ffDE 


n,„, a 


2 



"The abbreviations used in Table II and Fig. 10. 

'The free parameters in each model. Note that the additional parameter h appearing in the data fits is not considered as a model parameter but is included in k 
when calculating AIC and BIC. 



TABLE 11: Summary of the information criteria results 



Model 




AAIC 


ABIC 


A 


468.461 








w 


468.327 


1.866 


5.862 


aDE 


468.452 


I.99I 


5.987 


GCG 


468.461 


2 


5.996 


CPL 


467.663 


3.202 


11.195 


HDE 


470.513 


4.052 


8.048 


RDE 


493.772 


27.311 


31.308 


ADE 


503.039 


34.578 


34.578 


DGP 


530.443 


61.982 


61.982 



Notes: The cosmological constant model is preferred by both the AIC and the BIC. Thus, the AAIC and ABIC values for all other models in the table are 
measured with respect to these lowest values. The inodels are given in order of increasing AAIC. 



This model has the lowest values of AIC and BIC in all the models tested, so AAIC and ABIC are measured with respect to this 
model; see Table II. 

To make a comparison, we refer to Ref. [19] in which Komatsu et al. give the best-fit parameters: D.,,, = 0.273 and li = 
0.702, for the flat ACDM model from WMAP 7-year data combined with BAO and Ho data. We find that our results with SN 
Constitution data are consistent with this result. We also stress that in the joint data analysis we have used the chi square of SN 
data jp^^ that is /i-free, instead of;^'^^ that is /i-relevant. For making a clear comparison, we also perform a joint analysis with the 




0.38 0.40 



FIG. 1: The cosmological constant model: likelihood contours at 68.3% and 95.4% confidence levels in the Q.„, - h plane. The left panel 
shows that when we use the SN Constitution data with h not marginalized, there will be an inconsistency between the data of SN and other 
cosmological probes. The right panel tells us that the tension will disappear when we use the SN Constitution data with h marginalized. 
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FIG. 2: The constant w model; likelihood contours at 68.3% and 95.4% confidence levels in the Q„, - w and Q„, - h planes. 

/j-relevant SN chi square and in this way we find the best-fit parameters; Cl,„ = 0.361 and li - 0.644, wi\hxl,ii, = 495.55, 
in accordance with the resuhs of Ref. [41]. Such a big xl,j„ implies that this way does not seem to be correct. Figure 1 shows 
the probability contours in the Q„, - h plane for the flat ACDM model. The left panel tells us that if we use the h-iele\antxlf^, 
a great tension will be brought between the SN limit and the BAO+CMB limit, and the right panel shows that the tension will 
disappear when considering the h-freex^^. Actually, in Ref. [42], Wang and Mukherjee have argued that because of calibration 
uncertainties, SN data need to be marginalized over li if SN data are combined with data that are sensitive to the value of Our 
results further confirm this opinion. 

B. Dark energy models with equation of state parameterized 

For this class, we consider two models; the constant w parametrization and the Chevallier-Polarski-Linder (CPL) parametriza- 
tion. 

1. Constant w parametrization 
In this case, the equation-of-state parameter of dark energy is assumed to be a constant, so in a flat universe we have 



E{Z) = Vj^m(l + Zy + + + (1 - a„ - + 



This is a two-parameter model with the model parameters 0,„ and w. 
The best-fit parameters and the corresponding are; 



Q,„ = 0.272, w = -0.981, h = 0.703, xL, = 468.327. 



(32) 



(33) 



We plot the likelihood contours for this model in the Q,„ - w and 0,„ - li planes in Fig. 2. From this figure, we see that when 
the equation of state does not depend on redshifts, dark energy is consistent with a cosmological constant within Icr range. 
Comparing with the cosmological constant model, this model gives a ^ower xlu,,' ^lut due to one extra parameter it has, it is 
punished by the information criteria; AAIC = 1.866 and ABIC = 5.862. 



2. Chevallier-Polarski-Linder parametrization 



Now, we consider the commonly used CPL model [4], in which the equation of state of dark energy is parameterized as 

z 



w(z) - Wo + W„ 



1 +z' 



(34) 
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FIG. 3: The Chevallier-Polarski-Linder (CPL) model: likelihood contours at 68.3% and 95.4% confidence levels in the wo - Wa and 0„, - h 
planes. 



where wo and Wa are constants. The corresponding E{z) can be expressed as 



Eiz) = ^n,„il + z)3 + Q,(l + z)4 + (1 - Q,„ - Q,)(l + 2)3(i+u.„+H,) exp (-^j. (35) 

There are three independent model parameters in this model: = {0.,„, wq, Wa}- 

According to the joint data analysis, we find the best-fit parameters and the corresponding ;if^j|jj: 

Q„, = 0.265, wo = -0.847, w„ = -0.691, h = 0.716, xlin = 467.663. (36) 

We plot the likelihood contours for the CPL model in the wq - Wa and Q.,„ - h planes in Fig. 3. We note that the best-fit 
parameters of the ACDM model (Q,,, = 0.275, vvq = -1 and Wa = 0) still lie in the Icr regions of the CPL model, indicating that 
the ACDM model is fairly consistent with the current observational data. Since the CPL model has three free model parameters, 
it should have made considerable improvement in the fit, however, it gives a nearly equal contrasting to the constant w model 
(only smaller by 0.664). The differences in the information criteria with respect to the ACDM model are: AAIC - 3.202 and 
ABIC = 1 1 . 195. The information criteria result of CPL is worse than the w model (especially its ABIC value is very large). This 
implies that the CPL model is too complex to be necessary in explaining the current data, comparing with the simpler models 
such as the A model and the constant w model. 



C. Chaplygin gas models 

Chaplygin gas models describe a background fluid with p cc p " that is commonly viewed as arising from the d-brane theory. 
Moreover, these models may be able to unify dark energy and dark matter. We should have considered both the original (o- = 1) 
and the generalized Chaplygin gas models, however, the original Chaplygin gas model [12] has been proven to be inconsistent 
with the observational data [22], we thus only consider the generalized Chaplygin gas (GCG) model [13] in this paper. 

The GCG has an exotic equation of state: 

Pgcs = -7^' (37) 

rgcg 

where A is a positive constant. This leads to the energy density of the GCG: 

Pgcgia) = PgcgiO) Ja, + ^5^1 , (38) 

where A j = A/p^^t^(0). When A , - 0, the GCG behaves like a dust matter; when A j = 1, the GCG behaves like a cosmological 
constant. So, the GCG model is considered as a unification scheme of the cosmological constant and the CDM. In a flat universe. 
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FIG. 4: The generalized Chaplygin gas model: likelihood contours at 68.3% and 95.4% confidence levels in the - a and Aj - h planes. 



we have 



E(z) 



7^(1 + + a(l + zr + (1 - - a) {As + (1 - A,)(l + .)3(i+«))i/i+«. 



(39) 



This model has two independent model parameters: = {As, a}. The cosmological constant model is recovered for a - and 

= 1 - - A,(l -Q.r- ^b). 

The best-fit parameters and the corresponding are: 



A, = 0.758, = 0.003, = 0.701, ^^ 4^34^^ 



(40) 



We find that for the GCG model the xj„i„ value is the same as that of the ACDM model, which is an amazing coincidence. The 
best-fit value of a is so close to zero, implying that the ACDM limit of this model is favored. We plot the likelihood contours 
for the GCG model in the As - a and A , - li planes in Fig. 4. As a two-parameter model, the GCG performs well under the 
information criteria tests: AAIC - 2 and ABIC - 5.996. 

D. Holographic dark energy models 

Holographic dark energy models arise from the holographic principle of quantum gravity. The holographic principle deter- 
mines the range of validity for a local effective quantum field theory to be an accurate description of the world involving dark 
energy, by imposing a relationship between the ultraviolet (UV) and infrared (IR) cutoffs [43]. As a consequence, the vacuum 
energy becomes dynamical, and its density p^^ is inversely proportional to the square of the IR cutoff length scale L that is 
believed to be some horizon size of the universe, namely, p^ie ^ In this subsection, we consider three holographic dark 
energy models: the original holographic dark energy (HDE) model [5], the agegraphic dark energy (ADE) model [8], and the 
holographic Ricci dark energy (RDE) model [10]. We note here that, different from the previous several models, the holographic 
dark energy models do not involve the cosmological constant model as a subclass. 



1. Holographic dark energy model 



The HDE model chooses the future event horizon size as its IR cutoff scale, so the energy density of HDE reads 
where c is a constant, and Rei, is the size of the future event horizon of the universe, defined as 



r°° dt' r 

Reh- a \ — = fl 

Jt a J„ 



da' 



(41) 



(42) 
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FIG. 5: The holographic dark energy model: likelihood contours at 68.3% and 95.4% confidence levels in the Q.,„ - c and CI,,, - h planes. 



In this case, E{z) is given by Eq. (29), where the function Q.de{z) is determined by a differential equation: 



e(z) - 1 + 



l+z 



where a prime denotes d/dz, and 



3 

eiz) = - 



l + Ul+z)y 



1 + (1 +z)r 



(1 - Q.de(z)) + (1 + Wde(z))£ldeiz) 



with 



1 2 ; 



(43) 



(44) 



(45) 



The HDE model contains two independent model parameters: - {Q.,„, c). Solving Eq. (43) numerically and substituting the 
resultant ^deiz) into Eq. (29), the corresponding E(z) can be obtained. 

For this model, we obtain the best-fit parameters and the corresponding ;(f^j.^j: 



Q„, = 0.285, c = 0.742, li = 0.684, xL, = 470.513. 



(46) 



Our results are generally consistent with those derived in previous works [31, 44]. We plot the likelihood contours for the HDE 
model in the Q,„ - c and Q.,„ - h planes in Fig. 5. The HDE model performs fine under the information criteria tests, with the 
results AAIC - 4.052 and ABIC - 8.048. It should be stressed that the HDE model does not contain the ACDM model as a 
sub-model, whereas other two-parameter models which perform better than the HDE model, such as the aDE, constant w, and 
GCG models, all nest the cosmological constant A and tend to collapse to the ACDM model once being up against the current 
observational data. 



2. Agegraphic dark energy model 

The ADE model discussed in this paper is actually the new version of the ADE model [8] (sometimes called the new ADE 
model in the literature) which chooses the conformal age of the universe 



as the IR cutoff, so the energy density of ADE is 



f dt' r da' 

— = TT^ (47) 
Jo a Jo Ha'^ 



Pde = 3n^Ml,Ti-\ (48) 
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FIG. 6: The agegraphic dark energy model: likelihood contours at 68.3% and 95.4% confidence levels in the n - h plane. 



where n is a constant which plays the same role as c in the HDE model. 

As the same as the HDE model, E{z) is also given by Eq. (29), where the function D.de{z) is governed by the differential 
equation: 



1 + z \ n 



where the form of e{z) is the same as Eq. (44), in which 



Wde(z) = -1 



2(1+z)VQ^ 
3n 



(49) 



(50) 



Following Ref. [46], we choose the initial condition, ^deizmi) - n^i^ + ZiniT'^ l'^, at z,„, - 2000, and then Eq. (49) can be 
numerically solved. Substituting the resultant ildeiz) into Eq. (29), the function E{z) can be obtained. Note that in this model 
once n is given, by solving Eq. (49), Q,„ = 1 - Qi/c(0) - Q.,- can be derived accordingly. So, actually, the ADE model is a 
one-parameter model; the sole model parameter is n. 

For this model, we get the best-fit parameters and the corresponding ;^'^^|.^: 



2.755, h = 0.654, 



y2. 
A mm 



503.039. 



(51) 



This leads to Q.„, - 0.287. We plot the likelihood contours for the ADE model in the n - h plane in Fig. 6. As a single-parameter 
model, the ADE performs much worse than the ACDM model: its xl„j„ is greater than that of the ACDM model by about 30, 
and its AAIC = ABIC = 34.578. 



3. Ricci dark energy model 

In the RDE model, the IR cutoff length scale is given by the average radius of the Ricci scalar curvature I*??! '''2, so in this case 
we have pde ^ "R- In a flat universe, the Ricci scalar is !R = -6{H + 2H^), and as suggested in Ref. [10], the energy density of 
RDE reads 

Pde = 3aMl,{H + 2H^), (52) 
where a is a positive constant. From the Friedmann equation, we derive 

E- = ^,„oe-'' + Q,oe-^-'" + « (5 ^ + 2^' j , (53) 
where x = Ina. Solving this differential equation, we get the following form: 

E{z) = J|^(l + zy + ad + zY + (1 - a - |^)(1 + zr-^l (54) 
\ 2 - ff 2- a 
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FIG. 7: The Ricci dark energy model: likelihood contours at 68.3% and 95.4% confidence levels in the £1„, - a and fl„, - h planes. 



This is a two-parameter model, and its free model parameters are: 6 - {Q„,, a}. 
For the RDE model, the best-fit parameters and the corresponding are: 

= 0.371, a = 0.313, /i = 0.643, = 493.772 



(55) 



We plot the likelihood contours for the RDE model in the 0,„ - a and Q.,„ - h planes in Fig. 7. Like the ADE model, RDE also 
performs very bad: AAIC = 27.31 1 and ABIC = 31.308. This conclusion is consistent with the previous work [31, 45]. 



E. Dvali-Gabadadze-Porrati brane world and related models 



The DGP brane world model is a well-known example of the modification of general relativity for explaining the acceleration 
of the universe. In this subsection, we consider two models: the DGP model [15] and its phenomenological extension (namely, 
the a dark energy model) [17]. 



1. Dvali-Gabadadze-Porrati model 



The DGP model arises from the brane world theory in which gravity leaks out into the bulk at large scales, resulting in the 
possibility of an accelerated expansion of the universe. In this model, the Friedmann equation is modified as 

3Ml, {h- - = p,„(l + zf +Pril+ z)\ (56) 

where - {Ho(l - Q.„, - fl^)) ' is the crossover scale. In this model, E{z) is given by 

E{z) = ^Q,„(l+z)Ua-(l+z)4 + av + VQ^, (57) 

where Q.rc - ^li^^lH^) is a constant. The flat DGP model only contains one free model parameter, Q,„. 
For the DGP model, the best-fit parameters and the corresponding are: 

Q„, =0.295, /i = 0.632, = 530.443. (58) 

We plot the likelihood contours for the DGP model in the Q.„i - h plane in Fig. 8. We see that the DGP model, as a single- 
parameter model, is even worse than the ADE model under the observational tests. ^i^x^„j„ is greater than that of the ADE model 
by about 30, and it yields AAIC = ABIC = 61.982, also much larger than all other models we considered. So, the fitting result 
shows that the DGP model seems to be inconsistent with the current observational data (see also Ref. [22]). 

What should be mentioned is that the DGP model could perform much better when considering the systematic errors of the 
SN la data. For example, using the MLCS2k2 light-curve fitter for the SNe la data, the authors of Ref. [47] found that the 
DGP model performs better than the standard ACDM model. Currently the SNe la measurement errors are being dominated by 
systematic rather than statistical uncertainties, and for the sake of simplicity we would not discuss this problem in this paper. 
See Refs. [18, 41-49] for detailed discussions of this issue. 
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FIG. 8: The Dvali-Gabadadze-Porrati (DGP) model: likelihood contours at 68.3% and 95.4% confidence levels in the £1,„ - h plane. 

2. Phenomenological extension of DGP: a dark energy model 

Inspired by the DGP brane world model in which the Friedmann equation is modified, Dvali and Turner [17] proposed a 
phenomenological dark energy model, sometimes called the a dark energy model, which interpolates between the pure ACDM 
model and the DGP model with an additional parameter a. In this model, the Friedmann equation is modified as 

3Mli ^H^ - j = p„,(l + zf + p.(l + z)\ (59) 

where a is a phenomenological parameter, and = (1 - 0„, - Qr)'''*"^^^//^'. According to this Friedmann equation, E{z) is 
determined by the following equation: 

E{zf = n,„(l + + Q.r{l+ z)* + Eizra - Q„, - O,). (60) 

So, this model is a two-parameter model, with the independent model parameters 6 = {0,,,, o'). Note that a - 1 corresponds to 
the DGP model and a - corresponds to the cosmological constant model. 

Our joint analysis shows that for the aDE model the best-fit parameters and the corresponding are: 

Q,„ = 0.276, a = 0.030, h = 0.702, xL, = 468.452. (61) 

We plot the likelihood contours for the aDE model in the 0,„ - a and Q.,„ - h planes in Fig. 9. We notice that the best-fit value 
of a deviates one evidently, implying that the DGP model is incompatible with the current observational data. The cosmological 
constant limit, a = 0, is consistent with this model within Icr range. Moreover, we find that the aDE model gives tht x^„i„ 
smaller than the A model under our investigation, and its information criteria results, AAIC = 1.991 and ABIC = 5.987, also 
indicate that the aDE model fares the best, except for the ACDM and w models, under the current observational tests. 

V. DISCUSSION AND CONCLUSION 

In this work, we have considered nine popular dark-energy cosmological models and tested them against the latest cosmologi- 
cal data. This includes observational data of SNe la from the Constitution compilation, BAO from the SDSS, the CMB "WMAP 
distance priors" from the WMAP seven-year observations, and the measurement of Ho from the HST. We have used the "strong 
inflation prior" that imposes a flatness prior, and explored dark energy models in the context of such a flat universe assumption. 
To assess the various competing dark energy models and make a comparison, we have applied the information criteria, both the 
BIC and AIC, in this analysis. 

For each model, we have outlined the basic equations governing the evolution of the universe, calculated the best-fit values of 
its parameters, and found its AAIC and ABIC values. Table I summarizes all the models under consideration and the parameters 
that describe each model. We have plotted the likelihood contours of parameters for all the models. The fit and information 
criteria results have been summarized in Table II. Note that since the cosmological constant model has the lowest values of both 
AIC and BIC, the values of AAIC and ABIC are measured with respect to this model. 
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FIG. 9: The a dark energy model (phenomenological extension of the DGP model): likelihood contours at 68.3% and 95.4% confidence levels 
in the Q.„, - a and Q.„, - h planes. 




FIG. 10: Graphical representation of the results in Table II: the values of AAIC and ABIC for each model. The order of models from left to 
right is the same as the order in Table II, which is listed in order of increasing AAIC. 

From Table I, we see that according to the number of parameters the models can be divided into three classes: the one- 
parameter models, including A, ADE, and DGP models; the two-parameter models, including w, GCG, HDE, RDE, and aDE 
models; and the three-parameter model, namely, the CPL model. If we only com^&re. x^„m^ we find that the A model is not the 
best one, and there are four models, namely, o'DE, w, GCG, and CPL models, a little bit better than the A model according to 
this criterion. However, if the economical efficiency is considered, one will find that the A model is the best one to explain the 
current data, since both the AIC and BIC values it yields are the smallest. Although the CPL model can fit the current data well 
and has the lowest it is too complex (it has three free model parameters) to be necessary, yet. In the two-parameter models, 
the w model performs best in explaining the current data. 

We can also classify these models in another way that whether the model can reduce to the A model. Some models, such as 
w, CPL, GCG, and o-DE, can reduce to the A model, but the other ones, namely, HDE, ADE, RDE, and DGP, can not. From 
Table II, we see that the models nest the A all perform well, whereas the models that cannot reduce to the A perform illy, except 
for the HDE model. The HDE model has the ability in explaining the current data nearly as good as orDE, w, and GCG that nest 
the A. Also, we notice that although the o'DE, w, GCG, and CPL models can fit the data well, they actually all tend to collapse 
to the A model with their best-fit parameters. 

Out of all the candidate models we consider, it is obvious that the simplest A model remains the best one. Following it are 
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a series of models that give comparably good fits but have more free parameters. These include the oDE, w, GCG, and HDE 
models, which have two free parameters; and the CPL model, which has three free parameters. We have shown that the aDE, 
w, GCG, and CPL models can reduce to the A model and their best-fit parameters indeed do so (to within Icr ranges). The HDE 
model is the sole one that can give a good fit but does not nest A. The last three models in Table II, RDE, ADE, and DGP, are 
clearly disfavored. They have fewer parameters than models like CPL, but they score poorly because they are unable to provide a 
good fit to the data. They cannot reduce to the A model for any values of their parameters. We provide a graphical representation 
of the IC results in Fig. 10 which directly shows the scores (in the AIC and BIC tests) the models gain. 

In conclusion, given the current quality of the observational data, and with the assumption of a flat universe, information 
criteria indicate that the cosmological constant model is still the best one and there is no reason to prefer any more complex 
model. This conclusion is in accordance with the previous work by Davis et al. [22]. We look forward to seeing whether this 
conclusion could be changed by future more accurate data. 
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